Approved For Release STAT 



Approved For Release 


2009/08/1 7 : 

CIA-RDP88-00904R00010010C 


k 


Approved For Release 2009/08/17 : CIA-RDP88-00904R0001 001 00045-1 




Third United Nations 
International Conference 
on the Peaceful Uses 
of Atomic Energy 


A/C OH F . 2 8/ P/34 6 
UOSIi 

May 1^64 

Original s HU S3 I AW 


Confidential until official release during Conference 

EXTRACTION P ROCESSES AND THEIR MATHEMATICAL DESCRIPTION 

Rozen A.M. , Bezzubova A. I. , Elatomtsev B.V. ,Khcrkhorina L.P., 

Nemirovaky A.M. , Nikolotova Z.I., Phushlenkov M.F. ,Reshetko Yu.V. 

Shuvalov O.N. , Teterin E.G. , Vasiliev V.A., Yurkin V.G. 

In the preaent paper which generalizes and develops papers 1-15 ^ 
the quantitative regularities of the extraction chemistry and 
engineering are studied and interpreted. 

1 .Extraction Equilibria 

jL General. The extraction equilibrium regularities are discuss- 
ed using the extraction of uranyl and other actinide nitrates by 
neutral organophosphoric compounds(for general discussion see 4 ))as 
an example. These systems show the electrolyte-nonelectrolyte equili- 
brium when the chemical bond of an extractant with a compound being 
recovered is necessary to overcome the electrostatic interaction in 
an aqueous phase. However, the bond should be sufficiently weak to 
permit stripping(the interaction energy -c -Accordingly , to 

interpret the extraction power of solvents the theory of chemical 
bond is required while due to the chemical interaction weakness when 
describing the dependence of equilibria on the extraction conditions 
it is necessary to consider the contribution of the Van-der-Waals 
(up to 3kcal/mole)and especially electrostatic interactions to the 
chemical potential in terms of the solution theoiy(the electrolyte 
solution theory for the aqueous phase processes while the multi- 
component nonelectrolyte one for the organic phase)* 

The extract ion, not complicated by the hydrolysis or formation 
of anionic complexes is described by :Me n+ +nA"+qS+hH 2 0=MeA n qS.hH 2 0 
where A is an anion, S, an extractant. Accordingly, the distribution 
coefficient for unhydrated solvates(A=NO,)will be: 

* =y/*=K(No 5 ) n t ; +1 (s ) q r f/ r c (1.1) 

where K is a thermodynamic distribution constant ;y and x, concentra- 
tions of a compound in aqueous (a*p*)and organic (o.p.)phases ,the 
round brackets denote concentrations , fa are activity coeffici- 
ents of ions in a.p*,of a solvate in o.p.and of an extractant , the dot 
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motes the concentration coefficients: j*=l at x=o or y=0;^ s =1 atN =1 
The terms of the equation define s^an extractabtllity ;(N0,) 1 V* + = 
processes in a.p., salting outjS^the effect of dilution due to sol?a- 
tion(in case of the extraction of several compounds it describes the 
competition for a free extractant ) { J^-cfcC.iB the effect of the dilu- 
ent peculiarity (weak Interactions in o.p.). These effects are to be 
discussed consecutively. 

Ej rtractability depends on the extraction power of the solvent 
and the forces retaining a compound in a.p. According to 4 ^ »^) 

V/'o (o) aq (i« s ) ) 

where / (0 )is the Rvalue at a zero concentration, K s , the solvation 
constant in o.p. The values of ^(0)so far determined only for HC1 and 
HN0 3 are very lowty =3.10~ 5 and 2.1x10- 2 , respectively )i.e. ions are 
strongly bonded with water*' .Therefore for the effective extraction 
a strong solvation is needed in o.p. 

Extraction power of phosphine oxides is maximum and diminishes 

on substitution of alkyl groups R by electrophylic ones(Ro,Ph,Cl) 18 » 1? 
To obtain the quantitative regularities it is necessary to improve 

the criteria of extractability and structure. The effective constant 
K-K^ 6 Kd l£ and the e.lectronegativity(EN)of groups, X, were assumed 
to be such criteria '.To improve the EN scale 18 'the infrared spectra 
of a number of organophosphoric compounds R^R^PO were investigated 
and the following scales were obtained(fig.2a,b) : 

a)Xp=4.0, EX=X RI +X R2 +Z R5 =6+0.024(w po -1l70) ( 1 . 2 ) 

h)Xj,=3.9, 5fx=6.6+0.021 ( a) ro -1170) ( 1 . 3 ) 

Uroup H OR Ph QPh CH 2 C1 CCl^ ph OPH CH g Cl C H p ClCH7ccT x 

Scale Bg X, Spectral X eff.(from extraction) ^ 

18) 1.65 2.0 3.0 2.4 3.2 2.85 3.35 2.7 3.5 3.5 27 i 

1.2 1.85 2.0 2.9 2.3 3.1 2.6 3.2 2.65 3.4 3.1 2*6 4*3 

1.3 2,2 2.253.0 2.5 3.172.8 3.3 2.8 3.473.2 2J3 4I2 

Lg K(the free energy)drops approximately linearly with the increase 
of EN,the number of radicals, n QR and <O pQ when alkyl radicals are sub- 
diminish es^ . et ^ er eal ones (fig.l )» (i.e.when the negative 0-charge 
_ ’ lg K=A-B £ X =•■ A r B 1 n QR = A^BgWpQ (1.4) 

x)Particularly for HCl,due to which it is more poorely extracted by 
TBP(Khc 1~2. 10'‘ 4 ,Kmjo^=0. 2) , though the bond strength of both the 
acids with TBP is clo3e.K StHC1 «=7 f K s mro. « 10. 

xx)The experiments were carried out using CCI4 for which di£«1 .Thel/jPu 
extraction isotherms were measured for the K calculation reliability, 

^ H was determined from the lgK temperature dependence, 

- 2 . 
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where constants A and B depend on the nature of the compounds extrac- 
ted. A stronger exrtractability decrease when introducing phenyl or 
chlormethyl radicals is approximately described(flg.1b) x 'by eq(l.4). 
with the help of effective (increased )EN,X ef =X fl + A X ef =X a ^aX;where 
X fl is EN of the alkyl radical substituted, AX=X-X ,/Wsconst 2 
( e .S.Xc C i 3 -Xa=1,35;X®Ji 3 -Xa=2, 6 , /W ^,9). 

The influence of the entropy or steric factors is one of the 
reasons effecting a difference between X ef and X as EN is responsible 
only for the change of the binding energy but not for all the free 
energy. In fact, the deviation from the general regularity is less for 
AH as compared to AF(fig.2c ,d) ,i.e. -aH=A 3 ~B 3 £ X=A 4 ~B 4 00 pQ . 

With the hydrocarbon chain elongation n c ,EN is slightly decreas- 
ed only for the first members of the series while the steric diffi- 
culties are increased, which results in the weak maximum of K obser- 
ved at n c =6-7(fig.1d,n c -s£n c /3).To describe the effect of n c i.e.the 
influence of the entropy-steric factors, the dependence lgK=p-rq(C/0) U , 
determined for ethers, can be used as the EN change cannot explain the 
constant decrease and C/0=2n c> As for the highly branched radicals the 
effective hydrocarbon chain length may be l>n c ,the general equation 
of the extraction power is 

-AF°/RT=lgK=A-B2 X ef -.ql=A 3 -B 3 £ X - ql' (1.5) 

The similar dependencies are expected in extracting by amines 
and acidic extract ants (in this case B<0);X ef differing from X even 
more owing to polymerization^. 

Salting out Theory'* Abased on the electrolyte solution theory, 
allows the effect of the salting out agent to be predicted, when the 
affective diameter of the extracted ions, the hydration numbers of 
ions and of the salting out agent are known. JJhen computing it is 
convenient to express the activity coefficients by Earned' s equation, 
and in case of a linear concentration dependence of j' by Rozen's one, 

isrC^e^salt.^y^^e'O^+^^^salt. , j 1 -*) 

For TJOg and Pucf^nit rates lg^J^^O^-0, 46+0. 116 I,i.e£=<^ VI =0.1l6 
The approximate values of Harned's coefficients are: 


Cation 

H + 

Na + 

ml 

Mg ++ ,Be ++ ,Ca ++ 

Al +++ 

for U(VI) 

0 

+0.06 

0.08 

0.033 

0.054 

for Pu(VI) 

0 

- 

0.13 

0.053 

0.078 


zj Fig.lajb.cjI-fRCO-PO.IItRO^R PO,IIIROR 9 PO,IVR,PO(.C ftJ oC.):fig. 
1»(R0ToPh0 PO;2,(CH 2 CICH 2 0) 3 PO,3 r (RO) 2 PhPO,4-(RO) 2 cfil,PO(Ph-pEenyl) ; 
fig . 2 : DHPhP-dihexylphenylphosphonate , DHCLMP-dihexylehloromethyl- 
phosphonate* 
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Effect of Diluents. The solvate formation results in the 
extraction dependence on the extractant concentration S(<x~S q ) 
whichis not affected by the diluent nature; the latter determining 
the value of the diluent parameter = -The activity coeffi- 

cients in multicomponent systems taking part in dit may be express- 
ed and interpreted by the properties of binary rystems with the 
help of the solution theoxy. the activity coefficients in binary 
systems diluent-TEP and diluent -uranyl nitrate solvate, are approxi- 
mately described by the equation with two constants: f \ 

d, 2 — 5*, Vj_/ Vi , while the data on the three-component systems 
by the equation with data on the binary systems^at C=0.Acid and 
water taken into account the expression c^iC becomes most complicated 
and contains 20 constants .describing all the possible binary inter- 
actions ; accordingly di£ depends on the nature of the diluent and 
the compound being extracted(fig.3a,b)as well as on the acid concen- 
tration (fig. 3c J 33 . The dit values derived from the data on the acti- 
vity coefficients of binary systems and from the extraction data are 
in agreement .The dit values are similar for all the diluents 
except CHCl^ , although the non-ideality sign in the TBP-diluent solu- 
tion is different and $ ^varies by 100 times.lt is explained by the 
similarity of the properties of TBP and solvate (with the higher non- 
ideality in solvate systems)due to which the interactions TBP-diluent 
and solvate-diluent affecting c lit in an opposite way are partially 
compensated, b^ 1 — 2b 21 . As compared to CCl^dit is increased (fig. 3a) 
due to the predominance of the interaction with TBP(hexane)or with 
solvate ( CgHg ) . The positive non-ideality of the TBP and solvate solu- 
tions in saturated hydrocarbons shows that the interaction between 
similar molecules is stronger than that between different ones, which 
is typical of the Van-der— ffaals forces. The strong negative non-idea— 

x) infr( i -y*)f«Pz(0«-2 6 (6 13 -2- (ft a „)7 + fz t f}C^i2'Q(i-2'f l )-6 z3 ft/tfl; C, 

is ternary constant, nos of components :1 -diluent ,2,extractant ,3, 
solvate* if , the volume fraction; V, the molal volume. 

xx)For ■fehe extraction of metal ndcroamountsCy^-Ojwith diluted extrac- 
tant (yp-Ojone obtains simple eqns,the form or which depends on the 
selected standard state of the solvate: dit* — a pure solvate. 

Il-a diluted solution in dry TBP,III~a solution in TBP .containing 
water and acid„E.g.lgc/i£l=2b 21 ..b3i ,ig di£II=2b 2 i -*>31+132 ,i.e. the devia 
tion from ideality is due to the extractant -diluent ( d 2 -) ) .solvate- 

diluent (h,.) and solvate-extractant (b-p)interact ions. 

xxx)For clL . .C^and CHCl-from extraction di£ I /dit I rrr = k/k p - f 1,55; 
2 ? 35;0,O2;from f :1,55;2,42?0,013. (die^cef ^(0) /fo{Q) 4l).NosV 
curves in fig.4a,b,d:1 .5 CCl^ diluent,2,6-CgHg^3,7hexane j4,8 CHCl^. 

HA 
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7 111 thS TBP - CHC1 3 3ystem results from the hydrogen bond formation 
as evidenced by the spectroscopic data 21 » 14 >and the high value of 
the heat of mixing(fig. 4d ).The main reason for the moderate negative 
non-ideality in systems with CCl 4 and CgH^and in CHC1- -solvate) is 
an athermal effect,i.e.the solution entropy increase due to the 
difference in sizes of the component molecules (fig. 4a, b) .However, the 
negative non-ideality is not confined to the athermal effect as the 
entalpy of mixing is ne^ative.Thus ,in these systems the "molecular" 
negative non-ideality(r=^atherm ^1 )exists x Wch may be interpreted 
by the unstable compound formation 5 '.The uranyl nitrate solvate- 
TBP system is also characterized by the negative non-ideality (b « 
-0.6 )and a high value of the heat of mixing (hj -=-77!hnoie)which 2 may 
be interpreted by trisolvate U0 2 (NC> 5 ) 2 .3TBP formation with K =1-2. 

The interaction of diluents with solvates of Th and Pu(lV) nitrates 
is similar to that with uranyl nitrate solvate, -but for Th solvate, 
the positive non-ideality is stronger in the systems with saturated 
hydrocarbons which results in the formation of two organic phases. 

The solutions of Th solvate in CC1 4 and in CgHgare close to athermal 
ones ip 50cal/mole)i.e. there is even a small molecular positive 

non-ideality. When substituting ethereal groups by alkyl ones(TBP-T0PO 
series)the activity coefficients of CC1 4 ,C 2 H 4 Cl 2 and CHC1, and.hence, 
the extractant-diluent interactions are markedly changed only for 
TBPO ' 5 '(fig*4c) . 

Th e Extrac tion Is otherm Calculation. In case of the concurrent 
distribution of U,Pu and HNC^macroamounts one obtains 8 ^: 

y lf f U S 2 5 y Pu =f Pu s2 5 y H = ( f IH +2f 2H } S ( 1 . 6 ) 

S=2S o/( 1+f 1 H +f 2H )(l + ' /T+8F S 0 = s o /(1+Bu/2) (1.7) 

where y and x are concn.of compounds in organic and aqueous phases, 

S and S Q are free and initial extractant concentrations 

f D = V c U N0 3 ’ f Pu=*Pu x Pu N0 3' f lrtH x H N0 3 * f 2H = ^2H x H N °3 
N° 3 =x h+ 2x u+ 4x + z i i s lt i ;B u =0,00247 S <> /(0,747 + O232S o ) f fl_- 0 _/ 2 . 

For TBP k 0 J^ r ’»7.1.)05-5«4(x o .x H ^«p u ) jtiB- , i , 9>k2if0jO0o4 _J f . 

component concentrations are related to the solvent volumes, then 
B H =e u =:0 » • 95 . 1 0° * ^ 2 5 ^ X U +X H //5 +X Pu ^-k 1 H =0 . 1 74 : k 2H =0 . 0002 . 

X l Pr ? bably ’ there is a weak acceptor-donor interaction between 

°£ benzene and. uncompleted Cl 3d-orbitals with TBP. 
xx)Cqns.(1 .6)are obtained by the combined solution of U.Pu and H1T0 
aistrl Dirt ion eqns(1#1) # 3 

Tt4 
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Extraction Cascades 

Calculation of Cascades. If the apparatuses are not numerous 
it is convenient to calculate graphically (fig. 5n)the macrocomponent 
distribution over the extractor stages by the method of successive 
approximations^ ^ . The distribution of elements present in micro- 
amounts is calculated analytically , fig. 6. The equations are obtained 
by the combined solution of equations of equilibrium y^c^x^and 
material balance which for some 1-stage is: 

L i_1 x i _ 1 +v i+i Y i+ 1 =I i x i +V i y ± , (2.1) 

where L and V are volumes of a.p.and o.p. The decontamination 
factor K^,in the extraction section atcx^n ^.1 determined by the 
extraction conditions at the first stage does not depend on the 
cascade(column)length. Conversely, in the scrubbing section atot'h'^ 

^-1 ,K d depends on the number of the stages;if a '=const. 

K d ^casc.~(1/a 1 n 1 )»1/(a'n') N 
If c*'n' >1 ,the scrubbing is not effective: 

K d,scrub.~ 1+1 ^ /n ^ +1 (°V n 5(« 2 n ') + *** ~ 2 (2.2) 

e.g.for Ru at «/= 0.2 and n' .10 we obtain K , =*= 2. If the extrac- 
tion is accompanied by a chemical reaction(the reduction stripping) , 
the element concentration in the solution at the i-stage which 
would t)e ^ 3: ^ 0 = x i _ 1 + ii(y i+1 --y i )in the absence of the reaction, will be 
x.jsxlj’ e“ k i i, k^the rate constant^t^ the time of the contact. 
Assuming y i =o< i x 1 ,we obtain the equation describing the process in 
the cascade: 

X i ( ° < i + e kiti ) = x i- 1 + ny i+ i^i. 1 + « i+ i x i+ i (2*4) 

If the reduction is eff ective(e k i t i>lo)so that the Pu( IV) concentra- 
tion in a.p. at the i-stage inlet ~ 0,a simple relationship 
between the concentration of the elements not reduced and the 
initial one(y^ +1 )is obtained 

y 1® y N+1^1*2*J^N ; y i =y N+1^i+1^i+2 # * # V 
where A ± =1 +e k i t i/ effective stage reduction coefficient. 

The above calculation method based on the stagewise contact pattern 
describes the process in the column satisfactorily(fig.5c ) .If two 

x) The values of the rate constants may be estimated from data 22 \ 
k=2000 [Fe][Pu^ + ] , provided that Pu4+ = Pu(lV)/(1+ZB-iiroi)=: Pu(IV)/(1+S), 


where B .is constant of complexing with nitrate ions. It is interest- 
ing to J note that too low rate constants are obtained when using Q 
found by the extraction method3 c /;at Fe=0.04m/1 and xjp2m/l,S=360, 
T/^alc.^min.wb* le 10sec. The extraction method overesti- 

mates the 8-values. exp. 
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actinide elements present in macroamounts and when a great number 
of calculations are needed(the static characteristic determination) 
EC should be used. We employed the programmes of direct (Ural l)and 
stagewise (Minsk I) computing. In the latter case the equations of 
equilibrium( 1 . 7 ),of material balance(2 . 1 )and those of mixing are 
used. The programme provides for the minimization of the deviation 
of the calculated inlet concentrations from the fixed ones by 
affecting the assumed values of the waste concentrations x^. Direct 
computing lies in the combined solution of equations relating the 
material concentrations to the value at all the cascade stages 1 5) . 
The microcomponent concentrations are computed after the actinide 
and acid determinations using an additional programme. 

Static characteristics of the extractor processes are determi- 
ned by a number of variables (L q ,L' ,V o ,x£,xg,S,etc. ), their influence 
may be approximately described as an action of one generalized 
parameter J? ,the degree of approximation to the limiting conditions 
(the flow theoretically minimum one ratio) 1 ^ 

? ^p^min^ 1 ^ « 1 , (2-6) 

(where y,y p are maximum operating and equilibrium U-c one ent rat ions 
in o.p.,V and V min> are operating and minimum o.p. flows). It is seen 
from the illustrative calculation of the Purex extraction-scrubbing 
column(the impurity effect on the Zr and Pu microamounts was 
neglected) (N=4 ;N'=1 ,1=76,5 , l/=56,2, V=363, xg=1.8m/l, xg=0.95m/l, 

(:i H ) '=f/ 1 » T o =1 - 2l,1 / l) * The data on U 8) ,HNO^ a) ,Pu 5) ,Zr 7b) and Ru 23) were 

used .The deconta mina tion from Ru was estimated by the separate 

distribution of tri~ and dinit ro so -nit rates RuT and RuDCcX- =c m a +C 

Ru T T D D 

where C is a fraction of species in a.p$)ata 2 ^ were recalculated 
for U-solutions(fig.7g) ;it was assumed that c2=f (x^)=0.025,and in 
the scrubbing section C T =1 , 0 (^= 0 ^. The flows L ,V and concentrations 
x£,T were varied.As it follows from figs. 7 , 8 with the y growth the 
uranium zone is expanded, the accumulation of Pu and Ru increases 
(the accumulation decrease at^>99^ results from the species losses 
due to the insufficient stage number), the waste concentrations of U 
and Pu and K d from Zr and Ru I ^rise.Fig.8 also shows that the static 
characteristics of the extractor are defined by the approximation 
dsgree to the limiting conditions while the way of approximation 

x jAstt,, >o<p.RuT is extracted completely and Kg^i/cS ~ 50 does not 
depend onojhowever the RuT accumulation is very large(fig„7e)and if 
the reaction RuT-RuD in the presence of U is rapid enough ,K & 
increases with . iki Ra 
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is of little importance. This conclusion is also true of the 
extraction stage number change. 

Optimization of Conditions. A possible way of optimization 
is a maintenance of the maximum K^(high/?)at the permissible U,Pu 
losses. Introduction of the optimization criterion is an alternative 
method. Assuming the conditions providing for the maximum K d with 
the minimum loss of U and Pu and minimum Pu accumulation to be 
optimum, a criterion is 

0 pt.=^rKju/(xJ u ) r (q pu ) 3 (2.7) 

where m,n,r and s are coefficients characterizing the significance 
of each factor when processing various fuels*). In many cases all 
these factors are equal ,m=n=r=s=1 (fig. 9d, curve l).’,7hen processing 
natural U the Pu accumulation is of no importance f s=0(curve2). If 
the highly active material is processed the importance of the 
decontamination increases and m=n=2 may be assumed(curve 3)The 
range of the optimum conditions is y =0. 90-0. 99 ;the greater is the 
importance attached to the decontamination, the higher isy 0 p t The 
use of conditions close to limiting ones is possible only with a 
sensitive systems of control(e.g.with that retaining the uranium 
front location in the extractor). 

Extr act ion Dynamics. If dy^dx^zrconst , the nonstationary mass- 
transfer eq. may be solved either approximately by the similarity 

method^ ,2 ^or exactly by the operational one^^.At t^t 

x 0 -x(z,t)^(x 0 -x p )(1-e" t/,t o), ° (2.8) 

where Xp(z)=x(z,°°)is the equilibrium concentration, while the relaxa- 
tion time t Q (the time constant)is defined from the relation"' 


t =removed material/initial transport=M/j 


(2.9) 


where j=V y p (x 0 )=L x 0 an,H=M 0 -U p 4i^ ol (x o -x)+ff ® ol (y o -y) ; g , holdups. 

More accurately t =(1-A)t ,the correction A is estimated from the 
24) u 

accurate solution .As in case of any separating cascades the 

similiarity of the stationary and non- stationary concentration 

distribution corresponds to eq. (2. 8). The degree of approximation 

x) Introduction of the coefficients of decontamination from Ru and 
Zr as separate factors is useful owing to the different effect 
of acidity, the scrubbing temperature, etc .on them. 

XX) + n a8 !T v)M + K+T,)&L i -zyr y f^ ~-0, y P =«x ; E=z/h 

5=z/HTU ; H=N ; 77,11111 time(=HTU holdup/flow) ;x,y, averaged x,y. 
e * S. 7 < /*c/*7/H *X q /£ H) a p = ]/[U*£) e £ "ij *x 0 e SZ ; £ = uni. 

ILL ° - 8 ~ 
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to equilibrium ^=(x o -x)/(x o -ic p )=1~e'' t//t o depends on time rather 
than on the coordinate. Time of equilibrium setting in ±s(3-5)t Q , 
which corresponds to ^=0.95~O w 99.The starting period curve found 
experimentally for the columm of /^200mm,H=:2 ,8m, with Rashig's rings 
of 25x25x4, 5mm(x^=35g/l,x^=0.5ni/l,n=V/L=1 ,<x=1 )is close to the 
exponent (fig. 9a )while the relaxation time is close to that calculated 
from eq.(2.9) . 


L=230 1/hr; r coW =1 4 , 6min L=3C01/hr ^ column =1 1 ,2min 


Pulsation intensity, J=af, mm/min 
T5 4x78 4x1(58 0 


4x34 


4x108 


HETS, rn 1.11 
t 0 ,min,calC.4.35 
t Q , min, exp. 3. 3 
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When starting from a blank solution tne expression for the 
is close to the exponent(2.8)(t n =M/j n ,U=2° ol 'I+2 °°}*, j n =LxJ 


x) 


outlet 

with delay. In case of nonlinear extraction isotherms (uranyl nitrate 
in wide ranges of concentrations and acidities) the process equations 
were solved with the help of E.C .Simultaneously the dynamic charac- 
teristics of the pulsed column of /zf 200mm, 4 ,2m high with a packing 
of 15x15x1mm were experimentally determined while starting up from 
a" solution" (x£=250g/l,x£=1.55m/l, yJ +1 =y^=0.2g/l,T o =0.8m/l,L=758Ar, 
V=315 fi/hr)and. from an acid(x^=0, x^=2m/l)with I=700mm/mln.The start- 
ing period curves for uranyl nitrate are close to the exponents 
(fig. 9c, f )w±th the lag timeATwhich increases with moving away from 
the bottom of the col umn (it being the disturbance source, i.e. not 
saturated Q.p).3y contrast, in starting up from an "acid"when the 

concentration front moves down the column the timeA^a/7/ is increas- 

c 

ed from the top to the bottom(fig.9e,f).M p ~19kg is removed, j Q =27kg/ 
hr and eq 2.9 gives t o =0.7hr>t o er P *«0.4hr.When starting up from an 
acid Mp=7kg, j 0 =17,5 and t° alc *=0.4h.r(t® xp ^O.3hr). 

The results of the transition process computations are close 
to the experimental data 

Masstransf er 

The Driving force (D.f.)With non-linear extraction isotherms 
the boundary concentrations x^ and y^=<*i x £ carmo ' k excluded by 
a simple procedure from the expressions for the flux through the 

x) fa'idXi/dt + 8y ti dyi/dt = L(*i- t-Xi)+V(fcM -ytj) Qi ,i-stage holdup. 
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diffusion Boundary layers x \the standard masstransfer equation 
becomes incorrect. In case of the nonelectrolyte distribution the 
task may be solved using the irreversible process thermodynamics 
according to which the diffusion flux is proportional to the chemical 
potential gradient: j=-D 0 cd/'/dz=B(a-a i ), where 5=5^5^ « (Q/jf) (a verage^ , 
da/dz ~(a-aj)/d\ With the equilibrium condition a =ka x one obtains the 
general masstransfer equation with the D.f .expressed in terms of the 
activities *3=Ay(a*-a y )=A x (a x -ap ; VA y =1/B +K/S x ;a*=Ka x ; A X =KA ; 

An alternative useful method for the electrolyte extraction (when 
it is difficult to average B)is the calculation of the boundary 
concent rations, let y p =F(x) and x-x ± =AX, then y 1 =F(x i )=F(x-Ax) and 
the flux is j=B x< ax=flyy f£ x- ax )- y] .Knowing F and B one can obtain 
AX,X i ,y i and j(or directly j , substituting^ by j/fl x ) from this 
equation.lt is possible, e.^. to expand F(X-AX)in a series on 
A^=d/B x power.lt gives j=Ky (y p -y) ,1/K y =1/fl y+ (y' + The correction 

A~n(y p -y)/(yp + &j/&y) depends on B y , therefore with (y p -y) D.f. 
the resistances are not quite additive. 

Masstra nsfer in a real apparatus is complicated by longitudlal 
axial mixingCam) By the macroflows, and by transverse non- 
uniformity (tn). All these effects may be approximately described by 
a longitudinal diffusion with the effective coefficient D f =D (p +D +D. 
resulting in the HTU increase; 6 am n 

h= V D ef.Mf.M% ( i + \y sh k +h D* m) (3.13) 
where h^ is HTU defined by the masstransfer coefficient ;w the velocitj 

h k .^s, W» x , 

The transverse non-uniformity(inconslstency of the flow rates 
over the column cross-section)is of a particular importance as it is 
the principal reason of apparatus effectiveness lowering with a dia- 
meter and height incfease(h may be 5 and more times higher then h^. 

The diffusion correction h D =D/w is large with low rates, therefore 
the effectiveness of columns of a large diameter increases with the 
flow rates(fig.11a,h tn =h-h k ; h k =HETS of column 25mm in dia.If D is 
decreased with the increase of the hole numbers and D x drops with 
sectioning and in the latter case D y =D x =1 ,5m 2 /hr, then curve 1 cores- 

x) j =S j)( (x-Xi )=B y (y t -y) , where B=D/< f , the -.masstransfer coefficient ; 

diffusion boundary layer thicness p; , hence j=K( oqx-y) ;l/K=l/l l+VB*) 
D.f .is(a l x-y)>(y p -y), 

xx) more accurately where y^-d n F(x)/dx n 

xxx)more accurately b pic f i H*r)-l)h K /h] +h py [i +(«. n 1) hp X /h]/[J-(«n-I)h py /h] 

3 $ 6 
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ponds to Dy=5.5m 2 Ar > D x =3.5,the o.p.t.n.is greater than i.n. a.p.). 

The pattern peculiarity is an anisotropy, a large difference in the 
diffusion coefficients in the longitudinal (D „)and transverse(D = 

©I J. 

' D T + ®am' direction3 : ®ef ,- D i Dj - * As " the tranverse non-uniformity 
decreases with the mixing intensification one may expect 11b ^that 
^th^l/^x ^ where D+vi-fc increased with the column diameter, Cy^d 2 , 
or exp d c )and D ef ^x+C^/B^ "'so that 1 )D gf> and D^should come nearer 
when turbulized.(This assumption is confirmed by fig.11c).2)Oppo- 
site to the small columns, in large dia columns the pulsation may lower 
D ef. due to the transverse non-uniformity decrees e(D tn =C 1 /(D o +C 2 d pa((f J) ) S 1 n(l 

C 1 is high). In fact, for 500mm, dia column at Vf=lm/h without pulsation 

* p p 

D =2m /h,with I=500mm/min D =1m /h ,for 100mm column D ==0,5 and 
0,85m /h,resp^as D T =D T +cldp)lt is important that the transverse non- 
uniformity is a hydraulic phenomenon therefore its study and elimina- 
tion is possible using hydraulic scale-up without the experiments 
with masstransfer' .Constructive measures result from the theory based 
on the channeling pattern 2 ^ 

Hn'M-'W+V a V 1 ? 

where AL/L is the flow fraction evading the masstransfer,l is the 
by-pass length. These measures consist in the diminishing of L by 
setting phase distributors and the decrease of £ by the apparatus 
sectioning, providing for the transverse flow mixing ;a good mixing is 
obtained by means of "louver" rotary plates 2 —^* 

The distribution of the dispersed phase in the packed column 
70Cttm in diameter was studied by hydraulic Bcale-up.With the pulsa- 
tion intensity increase the feed cone section expands while the 
number of feed sources required for the uniform distribution is 
decreased (fig. 12b). Flow density levelling off with pulsation was 
also observed for the distributors of a hydraulic seal type(fig.12c,d) 
The uniformity degree of the continuous phase motion may be controll- 
ed by an impulse method^^ ;it being important to scale up the density 
gradient that increases the by-pass length and, correspondingly, HETS 
in the uranium zone(fig.5c,also 2 ®^). 

Mechanism of the extraction intensification is often associated 
with the surface renewal, i.e. with the increase of the masstransfer 

x)Thus,the effectiveness of 5m in dia distillation plates with a 
directed motion was increased by 505^^). 

346 
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coefficient .The analysis of the elemental components of the process 
and of the experimental data on the sieve-plate and packed pulse 
columns 'showed that the effectiveness rise is due to the increase 
of the phase contact surface , while the effective coefficient of mass- 
transfer is even decreased, K^l ,fig.15a/( K £ =hoA;K R =do/d;K fl =i?/S? , 

V S / S o and V k A 0 are coefficients of effect ivlty, dispersion hSldup 
surface and masstransfer variation index "o"denotes the value before 
intensification; K e= K r K K m= K 6 K B K g ). m the first approximation the 
ic and k variations are secondary and are defined by the drop diameter 
decrease .Then the knowledge of the dependence l)of the drop radius 
on the intensification variahle( revolutions n, pulsation Intensity I) 
and 2)of the terminal velocity and the masstransfer coefficient on 
the drop diameter d is sufficient to advance a quantitative theory. 

The dependences of d on the intensity of I/d-v 1 ” 1 * 5 /G°* 5 type fo und 
for various apparatuses 1 1 * 12 ^correspond to the theoretical ones 25) , 
i.e.the drop size is determined from the general regularities of a 
dispersion in a turbulent flow. The terminal dropvelocity dependence u 
on its diameter is seen from fig. 14a(generalised coordinates, 

, r=c//v^ *d/d maXl s=dp<r/ / a^d/d min - t = = S/R = 
~dmax/dmin = VP ). The flooding flow rate w^ proportional to the terminal 
drop velocity sho'ald be decreased when dispersed ;the dependence of 
w £ on 1 should reproduce fig. 14a in a mirror image which is really 

ohserved(fig.15b,c),w f ^w^/fl+0,265(cl ef ) 2 j 12 ^ *), 

The increase of hold-up with the intensity may be described by 
Pratt s eq.as the characteristic velocity(v o --w. E ~U)depends on Ijhence, 
the dependences of S? on I and on the flow rates are similar and are 
defined by the degree of approximation to the flooding conditions 
y=w/w f .The masstransfer coefficient is somewhat decreased on 
dispersioning which together with the longitudinal diffusion growth 
accounts for the dependence of on I 11 '.The effectivity increase of 
packed-pulsed columns in two systems at I<2m/min is approximately 
described by the eq. 1^=1 +b E lJ f 12) (fig. I5d), where b E =1 .7(if a col /d pa £ k 
<7, then b ~0.25d f!ol /d pack ). At very high I HTU rises due to the’ drop * 
coalescence (surface decrease )and to the influence of axial mixing 11 ^. 
To choose the optimum conditions criterion jzf=w/h schould be used 11 ’ 2 ^ 

x;in the packed pulsed colu mn s the filling becomes close packed 
resulting in wS. decrease. 

3 4 6 1 
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Fig#1 Fffect of structure on "the extraction power 
of neutral organophosphoric compounds. 

(a,c) effect of substituting R groups by R0 ones; (b)total 
correlation; I (ROJjPOjII (RO) 2 RPO;III R0R 2 P0;IV n 3 P0; 

1 ( R0 ) gPhOPO ; 2 ( C 2 H 4 C1 ) 5 P0 ; 3 ( R0 ) 2 PhPO : 4 (RO ) 2 CC1,P0 

•-our data;o-calculated from data 1 ^ ' ;~scale 18 ' :*»scale(1 , 3) ; 

(d) effect of the length of the hydrocarbon chain, calculated 
from the data:1- ; ^ ; (phosphonates) ;5 .6- 1 ^ (phosphates) , 

Fig. 2 Electronegativity and its effect on the 
extraction. 

(a)scnle(l ,2) , X F =4; X C1 =3.0; X Br =2.8; (b)scale(l .3) ; X p =3.9; 

’ ^Br = ^ . 9 ; (c)(d)dependence of the binding energy 
( ^ - flH kcal/mole)onrXandW po , cm -1 , A,B -scaled . 2) (1 .3) ; 
0-TBP;TOP;1-methylphosphates; 2-0D0P; 3-T0P0; 4-(PhO) 3 PO; 
5-dihexylphenylphosphonate; 6-dihexylchloromethylphosphonate. 

Fig. 3 Effect of diluents and HROj on the extraction 
of metals by TBP. 

(a,b) extraction of microamounts of uranyl nitrate from 
2M NaNO-j . 

Diluents ; 1 -CCl^ ; 2-kerosene ; 3-benzene ; 4-CHCl^ ; a! -distribution 
coefficient in dry TBP; — ideal solution. 

(c,d) Extraction from nitric acid solutions. Diluent : 
benzene (data by Bleskina 1. A. ):!-“( VI) ; 2-Np(Vl); 3-Pu(IV); 

4-Kp ( IV ) ; 5-Th(IV) (5^=1 ,5M) ; a^asi-ideal solution. 

Fig. 4 Characteristics of the interaction of diluents 
with extractant and uranyl nitrate solvate. 
Activity coefficients in the systems diluent-TBP(a) and di- 
luent-solvate (b); indexed lines (0, A , a , x)_ j* .without 
indexes - ^ 2 or of athermal solution. Nos of 

curves: 1,5 CCl^ diluent; 2,6 CgHg; 3,7 hexane, 4.6 CHC1 5 . 

(d) thermodynamic excess functions of the systems; 

(c) effect of substituting R groups by RO ones; 

o-TBP; O -DBBP; A -BDBP; X-TBPO; (e) description of ternary 

systems according to the characteristics of binary ones; 

/ -calculated/. 
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Fig. 5 Calculation of the macrocomponent distribution 
over the extractor stages. 

(a) HNO^; (b) uranyl nitrate, g/l ( calculation with the 

volume changes being neglected); (c) distribution of uranyl 
nitrate and acid over the column height ( n and o-experi- 
mental; k and • -calculated); (d) operating line of nonsta- 
tionary (1) and stationary (2) processes. 

Fig, 6 Calculation of the distribution of elements 
present in microamounts for the extractors 
with various flow sheets. 

0C i - distribution coefficient over the i- stage ; & i = ot^n 
extraction factor, IT- Nos of stages; A-cascade extraction 
factor; R- recirculation coefficien when scrubbing; K-de- 
contaminat ion factor. 

Fig. 7 Effects of degree of approximation to the 
limiting conditions on the distribution of 
compounds over the extractor stages. 

(a) U,g/l; (b) HN0j f m/l; (c) Pu; (d) Zr; (e,f) RuT and Ru 
(without disproportionation ; values of ^ :1-0, 76; 2-0,90; 

3-0 , 99 ? 4-0 , 99 ; 5-0 . 9999 ; 6-0 , 999999 ; N-Nos of stages; 

(g) distribution coefficients of RuT and RuD(Xg,g/l;X H ,m/l) . 

Fig. 8 Static characteristics of the extractor and 
their optimization. 

Waste U,g/l,(a) and Pu( ^concentrations vsij ;1 ,2,3 ,-varia- 
tions of T or V,L q and respectively; 4-curve 2 with 
one additional extraction stage ;(c) factor of decontami- 
nation from Zr and Ru. 

1 -reaction RuT-*RuD does not proceed; 2-reaction is rapid; 
3-RuT is not accumulated; — curve 1 at N=6; (d)optimizatioi 
criterion vs n (1,2 -curves 1.2 without the RuT disproporti* 
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Fig. 9 Dynamics of the pulsed column of f! 200mm 

(starting period curves for uranyl nitrates 
in different column sections)* 

Starting from operating (a,b,c,d) and blank (e,f) solutions; 
(a,b) system with =n=1 at I=4x78mm/min, X^g/l; t,hr;H,m. 

_ . . 4V Flg * 10 Transition processes in the extraction column. 

f ! l0n ° f the concentratlon of uranyl nitrate (a,b) and 
HN0 3 (c,d) over the column height when starting from initial 
(a,c) or blank (b,d) solutions. Starting period curves for 
nitric acid when starting from initial (e) and blank(f) solu- 
tions; Xy.g/l- x h -M; t,hr; H,m. 

Fig. 11 Characteristics of transverse nonuniformity 
of the extractor, 

(a) Load characteristics of the packed column of ^400mm 

( I- a distributor with 50 hole; Il-with 200holes;3-section- 

lnS) ?fi n ?ihf * 25m “ (<3Urve 4) 111 case of uranyl nitrate extrac- 
tion °’ nD '; h,HETS,m;V ,m/hr; 

(b) applicability of the diffusion pattern; 

(c) D e f(1) and Ap ( 2 ) f cm 2 /sec vs revolutions n min7 1 of the 
column of $6 200mnr^ * 

Fig. 12 Hydraulic scaling up of dispersed phase 
distribution in the extraction column. 

(a) Relative dispersed phase distribution density q/q(sir.gle 
feed source in the centre and near the wall of column; r, mm; 
I,mm/min); (b) effect of pulsation intensity I ,mm/min , on the 
necessary feed source number (1) and the feed cone diameter, 
dcm(2); (c) dispersed phase distribution over the column 
cross-section, 1 ,1=0} 2. 1=1 500. 

Fig. 13 Masstransfer in case of nonlinear extraction 
isotherms . 
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1 4 Hydrodynamics of drops. 

(a) Terminal velocity u of a drops vs its diameter (genera- 
lized coordinates, calculated from data 295 ). Systems: water+ 
A - tetrabromethane; a - dibromethane , a - ethyl-bromide; 

• - brombenzene, o- nitrobenzene, +nitrotolyene, a - chlor- 
benzene, (b) Drag coefficient vs Reynolds number; --- rigid 
spheres. 


Fig. 15 Effect of the pulsation intensity I on the 
extractor characteristics, 

K R» K it and ( column of / 100mm), The limiting loads of 

sieve plate (b) and packed (c) columns. 

(d)HETS; o,^ ,+ - columns 100,200 and 4O0mm in diameter, 
respectively; packing of 1 5x1 5x0.6mm 12 ’ 26) ; - correlation 125 
I ef =Ie c ( f - 10O >, m/mi*. I=af . 
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